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SUMMARY 


Constant  amplitude  tension-tension  fatigue  tests  were 
conducted  on  AS4/3501-6  graphite  epoxy  (02/e2/-B2)s  laminates, 
where  0  was  1 5,  20,  25,  or  30  degrees.  Fatigue  tests  were  conducted 
at  a  frequency  of  5  Hz  and  an  R-ratio  of  0.1.  Dye  penetrant  enhanced 
X-radiography  was  used  to  document  the  onset  of  matrix  cracking  in 
the  central  -6  degree  plies,  and  the  subsequent  onset  of  local 
delaminations  in  the  0/-0  interface  at  the  intersection  of  the  matrix 
cracks  and  the  free  edge,  as  a  function  of  the  number  of  fatigue 
cycles. 

Two  strain  energy  release  rate  solutions  for  local 
delamination  from  matrix  cracks  were  derived:  one  for  a  local 
delamination  growing  from  an  angle  ply  matrix  crack  with  a  uniform 
delamination  front  across  the  laminate  width,  and  one  for  a  local 
delamination  growing  from  an  angle  ply  matrix  crack  with  a 
triangular  shaped  delamination  area  that  extended  only  partially  into 
the  laminate  width  from  the  free  edge.  Plots  of  Gmax  vs.  N  were 
generated  to  assess  the  accuracy  of  these  G  solutions.  The  influence 
of  residual  thermal  and  moisture  stresses  on  G  were  also  quantified. 
However,  a  detailed  analysis  of  the  G  components  and  a  mixed-mode 
fatigue  failure  criterion  for  this  material  may  be  needed  to  predict 
the  fatigue  behavior  of  these  laminates. 
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INTRODUCTION 


In  part  I  of  this  study  [1],  the  damage  observed  in  quasi-static 
tension  tests  on  (O2/02/-02)s  AS4/3501-6  graphite/epoxy 

laminates,  where  6  was  15,  20,  25,  and  30  degrees,  was  documented. 
In  addition,  the  stresses  responsible  for  the  matrix  cracking  and 
associated  local  delamination  observed  were  calculated  using  finite 
element  analyses. 

In  part  II  of  this  study,  fatigue  tests  were  conducted  on 
(O2/02/-02)s  AS4/3501-6  graphite/epoxy  laminates,  where  0was 
15,  20,  25,  and  30  degrees.  Constant  amplitude,  tension-tension 
fatigue  tests  were  conducted  at  a  frequency  of  5  Hz  and  an  R-ratio 
of  0.1.  Dye  penetrant  enhanced  X-ray  radiography  was  used  to 
document  the  onset  of  matrix  cracking,  and  the  subsequent  onset  of 
local  delaminations  at  the  intersection  of  the  matrix  cracks  and  the 
free  edge,  as  a  function  of  the  number  of  fatigue  cycles.  Two  strain 
energy  release  rate,  G,  solutions  were  introduced  for  local 
delamination  from  matrix  cracks.  These  G  solutions  were  used  to 
correlate  the  onset  of  local  delamination  in  fatigue  between  layups 
with  different  off-axis  angle  plies. 

EXPERIMENTS 


Materials 

As  documented  in  part  I  [1],  panels  were  manufactured  in  an 
autoclave  at  NASA  Langley  Research  Center  from  Hercules 
AS4/3501-6  graphite  epoxy  prepreg  using  the  standard 
manufacturer's  procedure  with  a  maximum  temperature  of  177°C 
(350°F)  in  the  cure  cycle.  One  panel  was  produced  for  each  of  the 
(O2/02/-02)s  layups,  where  0  was  15,  20,  25,  or  30  degrees.  Panels 
305  mm  (12  in.)  square  were  manufactured  and  then  cut  into  20 
individual  coupons,  127  mm  (5.0  in.)  long  by  25  mm  (1.0  in.)  wide. 
Each  specimen  had  a  nominal  ply  thickness  of  0.124  mm  (0.0049  in.) 
and  a  fiber  volume  fraction  of  67%. 

Experimental  Procedures 

All  static  tests  were  conducted  using  the  procedure  outlined  in 
part  I  of  this  study  [1].  Fatigue  specimens  were  placed  in  the 
hydraulic  grips  of  a  hydraulic  load  frame  such  that  each  specimen 
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had  a  gage  length  between  the  grips  of  approximately  76  mm  (3.0 
inches).  A  25  mm  (one  inch)  gage  length  extensometer  was  mounted 
on  the  specimen  midway  between  the  grips  to  measure  nominal 
strain.  As  specimens  were  loaded  quasi-statically  to  the  mean  load, 
a  plot  of  load  versus  nominal  strain  was  generated  using  an  X-Y 
recorder.  The  extensometer  was  then  removed,  and  the  specimens 
were  loaded  at  a  cyclic  frequency  of  SHz  and  an  R-ratio  of  0.1.  Tests 
were  run  at  several  maximum  load  levels  below  the  critical  load 
required  to  initiate  damage  in  the  static  tests.  The  cyclic  loading 
was  stopped  at  fixed  cyclic  intervals  and  a  zinc  iodide  solution  that 
was  opaque  to  X-rays  was  placed  on  the  edge  using  a  hypodermic 
needle.  A  radiograph  was  taken  of  the  specimen  width,  with  the 
specimen  loaded  at  the  mean  load,  to  document  the  presence  of 
matrix  cracking  and  delamination. 

Table  1  lists  the  maximum  cyclic  loads  and  the  schedule  used 
to  obtain  radiographs  for  each  cyclic  load  level.  For  example,  tests 
run  at  maximum  cyclic  loads  of  5000  and  6000  pounds  were 
interrupted  every  1000  cycles  to  obtain  a  radiograph  of  the 
specimen.  Tests  run  at  a  maximum  cyclic  load  of  3500  lbs  were 
interrupted  every  1000  cycles  for  the  first  10,000  cycles;  every 
2,000  cycles  between  10,000  and  20,000  cycles:  every  5,000  cycles 
between  20,000  and  50,000  cycles;  and  every  10,000  cycles 
thereafter,  to  obtain  a  radiograph  of  the  specimen.  Loading  was 
continued  until  a  radiograph  clearly  indicated  the  presence  of  a  local 
delamination  associated  with  a  matrix  crack  near  the  free  edge. 
After  the  fatigue  specimens  were  removed  from  the  hydraulic  grips, 
the  edges  were  polished  and  examined  in  the  light  microscope. 
Photographs  of  the  polished  edges  were  taken  through  an  optical 
microscope  to  document  the  damage  observed  on  the  straight  edge. 

Initially,  some  coupons  were  coated  with  a  thin  film  of  water- 
based  typewriter  correction  fluid  on  either  edge.  After  a  few  fatigue 
tests  were  performed,  however,  it  was  apparent  that  the  onset  of 
matrix  cracking  and  local  delamination  was  difficult  to  detect 
visually  on  the  edge  with  this  brittle  coating.  Hence,  most  of  the 
fatigue  coupons  were  tested  without  coating  the  edges.  It  was  also 
easier  to  see  the  formation  of  the  matrix  cracks  in  the  radiographs 
if  there  was  no  coating  on  the  edge,  because  the  dye  penetrant  used 
to  reveal  the  damage  was  more  likely  to  accumulate  near  a  coated 
edge. 
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Experimental  Results 


Sialic..  Tejsts 

The  results  from  the  quasi-static  tests  were  documented  in 
part  I  of  this  study  [1].  The  radiographs  showed  triangular-shaped 
shaded  regions  bounded  by  the  free  edge  and  matrix  cracks, 
indicating  that  a  local  delamination  was  present.  Edge  micrographs 
showed  that  local  delaminations  were  present  in  both  e/-9 
interfaces,  extending  in  the  same  direction  from  curved  matrix 
cracks  in  the  central  -6  degree  plies. 

EaliQus.  Tssta 

The  same  type  of  damage  seen  in  the  radiographs  and  edge 
micrographs  of  specimens  taken  after  the  quasi-static  tests  was 
present  in  the  specimens  that  were  tested  in  fatigue  [1].  However, 
unlike  the  quasi-static  tests,  the  sequence  of  damage  was  evident  in 
the  series  of  radiographs  taken  during  the  cyclic  loading.  For 
example,  fig.  1  shows  radiographs  taken  immediately  after  the  onset 
of  matrix  cracking  (fig. 1a),  and  the  subsequent  onset  of  local 
delamination  from  the  matrix  crack  (fig. 1b),  for  a  (02/202/-202)s 
laminate  that  was  subjected  to  a  maximum  cyclic  load  of  15,570  N 
(3500  lbs).  Similar  radiographs  were  obtained  for  the  other  layups 
tested  in  fatigue.  At  the  number  of  cycles  corresponding  to  local 
delamination  onset,  all  of  the  specimens  had  at  least  one,  and  in 
most  cases,  several  matrix  cracks  in  the  central  -0  degree  plies.  A 
few  matrix  cracks  were  also  observed  in  the  +0  degree  plies  for 
some  laminates.  Table  2  lists  the  average  number  of  matrix  cracks 
in  the  gage  section  measured  from  the  radiographs  of  laminates 
taken  just  after  the  onset  of  local  delamination.  Local  delaminations 
were  present  in  the  0/-0  interfaces  near  some,  but  not  ail,  of  these 
matrix  cracks.  By  the  time  the  fatigue  tests  were  terminated, 
matrix  cracks  had  also  developed  in  the  +0  degree  plies  of  many 
laminates.  Delaminations  in  the  0/-0  interfaces  often  terminated  at 
these  +0  degree  matrix  ply  cracks  (fig.2). 

Figures  3-6  show  the  cumulative  number  of  cycles  to  the  onset 
of  matrix  cracking  and  subsequent  local  delamination  onset  for  each 
layup.  The  load  at  onset  of  damage  (matrix  cracking  and  local 
delamination)  under  a  monotonically  increasing  load,  determined  in 
part  I  of  this  study,  is  shown  on  the  ordinate.  Figures  7-10  show  the 


4 


cumulative  number  of  cycles  to  local  delamination  onset  minus  the 
number  of  cycles  to  onset  of  matrix  cracking.  Hence,  figures  7-10 
represent  the  number  of  cycles  to  delamination  onset  once  the 
matrix  crack  is  present. 

STRAIN  ENERGY  RELEASE  RATE  ANALYSIS 
Uniform  Local  Delamination 

In  reference  [2],  an  equation  was  derived  for  the  strain  energy 
release  rate  associated  with  a  uniform  local  delamination  growing 
from  a  90  degree  matrix  crack  (fig. 11a)  in  a  multidirectional 
laminate  using  the  classical  definition  for  G  as 

(1) 

2  dA  ' 

where  dC  is  the  incremental  change  in  laminate  compliance 
associated  with  an  incremental  change  in  surface  area,  dA,  created 
by  the  local  delamination.  This  solution  assumed  that  the  composite 
consisted  of  a  laminated  region  and  a  locally  delaminated  region  in 
series,  resulting  in  the  following  expression 

- 3-1  ,2) 

2mw2  ‘ttoELD  tELAM* 

where  P  is  the  applied  mechanical  load  on  the  laminate,  w  is  the 
laminate  width,  and  m  is  the  number  of  delaminations  growing  from 
the  matrix  crack.  The  parameter,  m,  would  have  a  value  of  2  if  the 
cracked  off-axis  ply  is  in  the  interior  of  the  laminate,  corresponding 
to  local  delamination  on  either  side  of  the  matrix  crack.  However,  m 
would  have  a  value  of  1  if  the  cracked  off-axis  ply  was  a  surface 
ply.  In  eq(2),  t  is  the  laminate  thickness,  Elam  is  the  laminate 
modulus  (as  calculated  from  laminated  plate  theory),  tto  is  the 
thickness  of  the  locally  delaminated  region  that  carries  the  load 
(i.e.,  the  laminate  thickness  minus  the  thickness  of  the  cracked  off- 
axis  plies),  and  Eld  is  the  modulus  of  the  locally  delaminated  region 
(as  calculated  from  laminated  plate  theory). 

As  shown  in  appendix  1,  eq(2)  also  applies  to  the  case  of  a 
uniform  delamination  growing  from  an  angle  ply  matrix  crack  as 
shown  in  fig. (11b).  For  the  laminates  in  this  study,  t  *  Nh  and  tLD  * 
(N-n)h,  where  N  is  the  number  of  plies  in  the  laminate,  n  is  the 


5 


number  of  cracked  off-axis  plies,  and  h  is  the  ply  thickness. 
Therefore,  noting  that  Nx  =  P/w,  eq(2)  may  be  rewritten  as 


G  =  JS_(_-! - 1-1  ,3, 

2m h  \(N-n)ELD  NElam) 

The  twelve-ply  (O2/02/-02)s  laminates  tested  in  this  study  may  be 
analyzed  as  six-ply  laminates  with  twice  the  measured  ply 
thickness.  Hence,  N=6,  n=2,  and  N-n*4.  Then,  noting  that  m=2  for  an 
interior  -6  matrix  ply  crack,  a  normalized  strain  energy  release  rate 
may  be  calculated  as 


=  _ i_) 

^2  4  UEid  6Elam/ 


(4) 


Bendina’Extension  Coupling  Influence  on  G 

The  value  of  Eld  in  eq.(4)  depends  on  whether  or  not  the  locally 
delaminated  region  consists  of  symmetric  or  asymmetric 
sublaminates.  For  example,  the  locally  delaminated  region  may  be 
thought  of  as  a  (0/6 )s  symmetric  subiaminate.  In  this  case,  the 
modulus  of  the  locally  delaminated  region  is  simply 

Eid  = - - -  (5) 

(N-n)h{Aii)-’ 

where  (An)-^  is  the  first  term  of  the  inverse  extensional  stiffness 
matrix  for  the  (O/0)s  sublaminate.  However,  as  noted  in  ref.  [3],  if 
the  locally  delaminated  region  consists  of  two  2-ply  asymmetric 
sublaminates,  then  the  presence  of  bending-extension  coupling 
should  be  reflected  in  the  modulus  calculation.  This  may  be 
accomplished  using  laminated  plate  theory  by  prescribing  an  applied 
mechanical  strain,  ex,  setting  the  applied  Ny,  Nxy,  and  My  equal  to  0, 
and  setting  kx  and  Kxy  equal  to  0.  The  Nx  that  results  from  this 
analysis  is  derived  in  appendix  2  and  may  be  used  to  determine  Eld 
as 
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(6) 


Eld  = 


2 

(N-n)h 


where  Di,  D2,  D3,  and  D  are  functions  of  various  terms  in  the  A,  B, 
and  D  matrices  as  defined  in  appendix  2. 

Any  difference  between  equations  5  and  6  reflects  the 
constraint  imposed  on  the  2-piy  asymmetric  sublaminate  by  the 
undamaged  laminate  and  the  grips  (if  the  sublaminate  extends  into 
the  grips).  Because  of  this  constraint,  the  sublaminate  is  not  able  to 
develop  a  kx  or  Kxy  curvature.  However,  under  tension  loading,  a  non¬ 
zero  Ky  curvature  may  develop  away  from  the  grips  in  the 
sublaminate.  For  the  type  of  sublaminates  that  form  as  a  result  of 
edge  delamination  [3-5],  this  influence  of  bending-extension 
coupling  may  be  significant.  However,  for  the  (O/0)t  sublaminates 
that  constitute  the  locally  delaminated  region  in  the  (O/0/-0)s 
laminates  with  -0  matrix  cracks,  the  influence  of  bending-extension 
coupling  is  small.  This  small  influence  is  shown  in  fig. 12,  where  the 
normalized  G  of  eq(4)  is  calculated  as  a  function  of  0  for  (O/0/-0)s 
laminates  with  -0  matrix  cracks  using  material  properties  in  Table 
3.  The  normalized  G  values  determined  by  calculating  Eld  from  eq(5) 
assuming  a  symmetric  (O/0)s  sublaminate  are  nearly  identical  to  the 
normalized  G  values  determined  by  calculating  Eld  from  eq(6)  using 
the  loading  prescribed  above  in  the  laminated  plate  theory  analysis. 
However,  unlike  many  of  the  sublaminates  formed  by  edge 
delamination,  the  (O/0)t  sublaminates  that  constitute  the  locally 
delaminated  region  in  the  (O/0/-0)s  laminates  with  -0  matrix  cracks 
will  also  exhibit  shear-extension  coupling.  This  shear-extension 
coupling  may  have  more  influence  on  the  calculated  Eld  than  the 
bending  extension  coupling. 

Shear-Extension  Coupling  Influence  on  G 


The  influence  of  shear-extension  coupling  on  the  calculated 
Eld  for  the  2-ply  asymmetric  sublaminate  may  be  determined  using 
laminated  plate  theory  by  prescribing  the  loading  outlined  above,  but 
setting  yxy  =  0  instead  of  prescribing  Nxy  =  0.  The  Nx  that  results 
from  this  analysis  is  derived  in  appendix  2  and  may  be  used  to 
determine  Eld  as 


Eld  = 


(N-n)h 


A  (A12D22  -  612622!  A  ,  (A12B22  -  A22Bi2\p 
\A22D22  -  B22B22*  A22D22  -  B22B22' 


(7) 
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This  Eld  value  reflects  the  shear  constraint  imposed  on  the  locally 
delaminated  region  by  the  undamaged  laminated  region.  The  shear 
constraint  results  in  a  greater  apparent  stiffness,  Eld,  and  hence,  a 
correspondingly  lower  strain  energy  release  rate.  Figure  12  shows 
the  normalized  G  of  eq{4),  calculated  as  a  function  of  0  for  (0/e/-e)s 
laminates  with  -0  matrix  cracks,  where  the  normalized  G  values 
were  determined  by  calculating  Eld  using  eq(7).  This  loading 
assumes  that  the  locally  delaminated  region  consists  of  two  2-ply 
(O/0)t  asymmetric  sublaminates  with  bending-extension  and  shear- 
extension  coupling.  As  shown  in  fig.  12,  the  normalized  G  values 
determined  by  calculating  Eld  using  eq(7)  are  significantly  lower 
than  the  normalized  G  values  determined  by  calculating  Eld  using 
eqs(5&6)  for  0  between  15  and  30  degrees. 

Uniform  Local  Delamination  vs.  Edae  Delamination 

As  noted  earlier,  the  delaminations  that  form  in  the  0/-0 
interfaces  of  (O/0/-0)s  laminates  are  localized  and  bounded  by  -0 
degree  ply  matrix  cracks.  Hence,  the  physical  evidence  dictates  that 
these  delaminations  should  be  analyzed  as  local  delaminations. 
However,  it  may  be  enlightening  to  calculate  the  strain  energy 
release  rate,  G,  for  edge  delamination  in  the  0/-0  interfaces  of  the 
(O/0/-0)s  laminates  and  compare  that  to  the  corresponding  G  for 
local  delamination.  The  strain  energy  release  rate,  G,  for  edge 
delamination  is  given  by  [4] 


G  ==-^(Elam  -  E  )  (8) 

where  e  is  the  nominal  applied  strain,  t  is  the  laminate  thickness, 
Elam  is  the  laminate  modulus,  and  E*  is  the  modulus  of  the  laminate 
once  it  is  totally  delaminated,  as  determined  by  the  rule  of  mixtures 


Eiti 

E‘=-i^J -  (9) 

t 

where  M  is  the  number  of  sublaminates  formed  by  the  edge 
delamination  and  Ej  and  tj  are  the  modulus  and  thickness  of  the  i 
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sublaminates,  respectively.  For  edge  delamination  in  the  0/-0 
interfaces  of  the  (O/0/-0)s  laminates,  eq(9)  becomes 

g.  ^  (10) 

3 

The  strain  energy  release  rate  calculated  in  eq(8)  will  depend  on  the 
assumptions  used  to  determine  the  (0/0)  and  (-0)  sublaminate  moduli 
in  eq(10),  which  determines  the  delaminated  modulus,  E‘,  in  eq(9). 

To  compare  the  different  possibilities,  eq(8)  was  rearranged 
to  yield  a  normalized  G  as 


-S- «  1|Elam  -  E'  ) 
e^t  2 


Figure  13  shows  the  normalized  strain  energy  release  rate,  G/e^t, 
for  edge  delamination  in  the  0/-0  interfaces  of  the  (O/0/-0)s 
laminates  as  a  function  of  0  where  the  sublaminate  moduli  were 
determined  assuming  they  were  either  symmetric,  asymmetric  with 
bending-extension  coupling  only  (B/E  Coupling),  or  asymmetric  with 
shear-extension  coupling  in  addition  to  bending-extension  coupling 
(B/E  &  S/E  Coupling).  Unlike  the  comparison  for  local  delamination, 
there  is  a  significant  difference  in  the  normalized  G  for  edge 

delamination  between  the  symmetric  sublaminate  case  and  the 

bending-extension  coupling  case.  This  difference  results  from  the 
effect  of  bending-extension  coupling  on  the  modulus  of  the  -0 
sublaminate  in  eq(10)  which  does  not  contribute  to  G  for  local 

delamination.  Furthermore,  when  the  influence  of  shear-extension 
coupling  is  included,  the  normalized  G  for  edge  delamination 

becomes  very  small  for  all  values  of  0. 

In  order  to  compare  normalized  G  values  for  both  local  and 
edge  delamination,  Hooke's  law  in  the  form  Nx  =  NheELAM  .  where  N*6, 
was  substituted  into  eq(4)  to  yield 


G  .  3Elam  j_j _ 1  I 

e2t  2  UEld  6Elam/ 


(12) 


Figure  14  shows  G  normalized  by  e^t  for  edge  delamination  and 
uniform  local  delamination,  as  calculated  by  eqs.(11)  and  (12) 
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respectively,  as  a  function  of  0  for  delamination  in  the  0/-0 
interfaces  of  the  (O/0/-0)s  laminates.  The  sublaminate  moduli  were 
determined  assuming  they  were  asymmetric,  with  shear-extension 
coupling  in  addition  to  bending-extension  coupling.  Figure  14 
indicates  that,  for  the  same  applied  strain,  the  strain  energy  release 
rate  associated  with  a  uniform  local  delamination  from  a  matrix 
crack  exceeds  the  strain  energy  release  rate  associated  with  edge 
delamination  for  all  values  of  0. 

Part'al  Local  Delamination 

The  local  delamination  growing  from  an  angle  ply  matrix  crack 
modeled  in  fig.llb  has  a  uniform  delamination  front  across  the 
laminate  width.  However,  the  local  delaminations  observed  in  the 
experiments  have  triangular  shaped  delamination  areas  that  extend 
only  partially  into  the  laminate  width  from  the  free  edge  [1]. 
Therefore,  an  analysis  was  developed  to  account  for  the  fact  that 
the  local  delaminations  in  the  (O/0/-0)s  laminates  did  not  extend 
across  the  entire  laminate  width,  but  instead  were  limited  to  the 
extent  of  the  matrix  crack  in  the  -0  degree  ply.  The  model  was  based 
on  the  assumption  that  the  regions  near  the  free  edges  containing 
the  partial  local  delaminations  act  as  a  system  of  sublaminates  in 
series,  just  as  was  assumed  for  the  uniform  local  delamination,  but 
this  series  system  was  assumed  to  be  in  parallel  with  the 
undamaged  laminate  in  the  interior  of  the  laminate  width  (fig. 11c). 

Figure  15  shows  three  possible  ways  of  visualizing  the 
delamination  growing  from  the  matrix  crack  in  the  -0  degree  ply.  Fig. 
15a  shows  a  triangular  shaped  delamination  that  intersects  the  free 
edge  at  90  degrees.  Fig.  15b  shows  a  triangular  shaped  delamination 
that  intersects  the  matrix  crack  at  90  degrees.  Finally,  fig.  1 5c 
shows  a  triangular  shaped  delamination  that  intersects  the  free 
edge  at  an  angle  of  20.  Eventually,  if  the  local  delamination  in  the 
0/-0  interface  reaches  a  +0  degree  matrix  crack  and  arrests,  the 
local  delamination  will  appear  to  be  an  isosceles  triangle  with  +0 
and  -0  angles  intersecting  the  straight  free  edge.  The  relationship 
between  the  matrix  crack  length,  d,  and  the  delamination  length 
along  the  free  edge,  a,  is  different  for  the  three  cases  shown  in  fig. 
15.  These  relationships  are  given  in  Table  4.  Hence,  the  shape  of  the 
delamination  front  must  be  known  in  order  to  relate  the 
delamination  length  on  the  free  edge  to  the  length  of  the  matrix 
crack. 
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For  a  local  delamination  from  a  matrix  crack,  the  compliance, 
S,  where  S  =  1/E,  for  a  laminate  with  a  given  local  delamination 
area,  A,  in  an  interface  whose  total  area  is  A\  may  be  represented 
as  (appendix  1) 

S -  Slam  _  A 

S*-Slam  a*  (13) 


where  S*  =  1/E*ld  ,  and  E*ld  =  (fLoELoVt  as  defined  in  ref. 2. 
Assuming  a  partial  local  delamination  is  present  at  a  matrix  crack 
on  both  free  edges  (fig. 11c),  and  noting  that  a’=ka,  where  k  is  given 
in  Table  4  for  the  three  geometries  shown  in  figures  15  a-c,  then 
A  =  maa'  =  mka2  and  A*  =  2mLka.  Substituting  for  A  and  A*  in  eq(13) 
for  any  of  these  three  geometries  yields 


S-1/-  1  - 

^IldEld 


1 

Elam 


(14) 


Recalling  from  ref .4  that,  for  a  system  in  parallel, 

2  dA 


(15) 


and  noting  that  for  the  partial  local  delaminations  shown  in  fig. 11c 
A  =  mka2,  where  k  is  given  in  Table  4  for  the  three  geometries 
shown  in  figures  15  a-c,  then  dA  =  2mkada.  Substituting  V  =.  wLt  and 
dA  into  eq(15)  yields 


G  =  .lMei(L)dE 
2  2mk'a1da 


(16) 


The  modulus  of  a  laminate  containing  a  delaminated  region  near  the 
free  edges  loaded  in  parallel  with  an  undamaged  region  in  the 
interior  of  the  laminate  width  is  given  by  [4]  as 

E  =  (e*  -  Elam)-^  +  Elam 

(a)  (17) 


For  the  partial  delamination  case  shown  in  fig.  11c,  A’  =  2mLa'  and 
(A')*  =  mwL.  Substituting  these  expressions  in  eq(17),  and  noting 
that  a'=ka,  where  k  is  given  in  Table  4  for  the  three  geometries  in 
fig.  15a-c,  yields 


E=(E*-ELAM)^a  +  ELAM 


(18) 


Unlike  E*  calculated  for  edge  delamination  in  eq(9),  E*  in  eq(18)  is  a 
function  of  the  local  delamination  length  along  the  free  edge,  a. 
Differentiating  eq(18)  with  respect  to  a,  and  substituting  into 
eq(16)  yields 


(19) 


Substituting  E*  =  1/S,  where  S  is  given  by  eq(14),  into  eq(19)  and 
differentiating  yields 


G  = 


_ 1 _ 

(Ba  +  C) 


__Ea _ 

(Ba  +  C)2. 


(20) 


where 

B  =  i(-J _ l-K-!-) 

2\tLDElD  iElam/'L'  (21) 


and 

'AM  (22) 


Although  the  relationship  between  the  matrix  crack  length,  d,  and 
the  delamination  length  along  the  free  edge,  a,  is  different  for  the 
three  cases  shown  in  fig. 15,  eq.(20)  is  independent  of  k,  and  hence, 
G(a/L)  is  independent  of  the  local  delamination  geometry  assumed. 
However,  as  noted  previously  for  eq(2),  G  in  eq(20)  will  depend  on 
the  assumptions  made  when  determining  Eld  in  eq(21).  Therefore,  B 
in  eq(21),  and  hence  G  in  eq(20),  will  depend  on  whether  eqs.  5,6,  or 
7  are  used  to  calculate  Eld-  f^or  the  remainder  of  this  paper,  it  will 
be  assumed  that  the  sublaminates  that  constitute  the  locally 
delaminated  region  experience  both  bending-extension  and  shear- 
extension  coupling.  Hence,  eq(7)  will  be  used  to  calculate  Eld- 
Eq.(20)  may  be  rearranged  to  yield 

g  =  £?l/ _ D _ _ e _ \ 


1  2 


(23) 


where  D  =  BL.  When  a/L  =  0,  then  eq.(20)  becomes 


G  = 


£^(2D\ 
2m  \c2j 


(24) 


Substituting  for  D  and  C,  normalizing  by  e^t,  and  noting  that  tLD=4h 
and  t=6h,  eq(24)  becomes 


J£.  =  M _ 1 _ \ 

e2t  2  UEld  6Elam*  (25) 

which  is  identical  to  eq(12)  for  the  uniform  local  delamination. 

Figs.16-19  compare  normalized  G  values  calculated  using 
eq(23)  for  partial  local  delamination,  as  a  function  of  delamination 
length  along  the  free  edge  normalized  by  the  laminate  length,  a/L, 
to  normalized  G  values  calculated  from  eq(12)  for  uniform  local 
delamination  and  from  eq(11)  for  edge  delamination  modeled  in  the 
0/-0  interfaces  of  (O/0/-0)s  laminates,  where  0  =  15,  20,  25,  and  30 
degrees,  respectively.  The  uniform  local  delamination  and  edge 
delamination  solutions  are  independent  of  delamination  length,  but 
for  the  partial  local  delamination  solution,  G  decreases  with 
increasing  delamination  length  along  the  free  edge.  Both  local 
delamination  solutions  yield  greater  normalized  G  values  than  the 
edge  delamination  solution  for  all  delamination  lengths. 

Flg.20  shows  G  normalized  by  e^t,  as  a  function  of  a/L  up  to 
a/L»0.1,  calculated  using  eq(20)  for  partial  local  delaminations  in 
the  0/-0  interfaces  of  (O/0/-0)s  laminates,  where  0  =  15.  20,  25,  and 
30  degrees.  Normalized  G  values  were  only  plotted  for  a/L  values  up 
to  0.1  because  the  magnitudes  for  small  matrix  crack  lengths,  and 
hence  small  delamination  lengths,  were  of  primary  interest  for 
comparing  to  experimental  data.  In  each  case,  G  decreases  slightly 
with  delamination  length,  unlike  the  solutions  for  uniform  local 
delamination  which  were  independent  of  delamination  length.  If  the 
partial  delamination  geometry  can  be  anticipated  and  linked  to  the 
matrix  crack  length  as  shown  in  fig.15a-c,  then  both  the  laminate 
toughness  as  a  function  of  fatigue  cycles  and  a  critical  matrix  crack 
length  would  be  needed  to  predict  the  onset  of  a  partial  local 
delamination  in  fatigue.  However,  as  shown  in  fig. 20,  if  the  critical 
matrix  crack  length  is  sufficiently  small,  and  hence,  the 
corresponding  delamination  length  along  the  free  edge  is  small,  the 
difference  between  the  uniform  and  partial  delamination  solutions 
will  be  insignificant.  In  this  case,  G  may  be  calculated  from  eq(2)  to 
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compare  to  experimental  data,  even  though  the  local  delamination  is 
not  uniform  across  the  laminate  width. 

Thermal/Moisture  Contribution  to  G 

The  strain  energy  release  rate  for  local  delamination  will  be 
influenced  by  the  presence  of  residual  thermal  and  moisture 
stresses  in  the  laminate.  A  uniform  local  delamination  G  solution 
that  includes  the  effect  of  the  mechanical  load,  the  residual  thermal 
stresses  that  exist  in  the  laminate  after  it  cools  down  from  the 
curing  temperature,  and  the  hygroscopic  stresses  that  exist  in  the 
laminate  due  to  the  presence  of  absorbed  moisture,  was  recently 
developed  [6].  Figs. 2 1-24  show  the  contribution  of  mechanical  load 
(M),  temperature  (T),  and  hygroscopic  (H)  moisture  weight  gain  to  G 
for  a  uniform  local  delamination  in  the  e/-e  interfaces  of  (O/e/-0)s 
laminates  (where  0  =  15,  20,  25,  and  30  degrees,  respectively) 
subjected  to  a  mechanical  load  of  3000  Ib/in.  and  a  AT  of  -156®C 
(-280°F).  As  shown  in  figs. 21 -24,  when  the  contribution  of  the 
residual  thermal  stress  is  included,  will  increase  above  the  G^^ 

value  obtained  from  equation  (2)  for  an  applied  mechanical  load 
alone.  However,  when  the  presence  of  absorbed  moisture  is  included, 
as  may  occur  from  exposure  to  the  ambient  air  (appendix  3),  the 
thermal  stresses  will  relax  and  reduce  the  strain  energy  release 
rate,  G^+T+h^  with  increasing  AH.  Similar  behavior  was  observed  for 
edge  delamination  [7].  However,  the  magnitude  of  the  thermal  and 
moisture  effects  is  much  more  significant  for  edge  delamination. 
The  laminates  tested  in  this  study  had  a  moisture  weight  gain,  AH, 
between  0.2  and  0.3  %  (appendix  3).  For  the  (O/0/-0)s  laminates 
shown  in  figs.  21-24,  this  corresponds  to  a  5. 8-7.7%  increase  in 
qM+t+h  compared  to  the  G^  solution  for  mechanical  loading  only. 
The  difference  in  G^+t+h  and  G*^  is  even  less  at  higher  mechanical 
load  levels  (Table  5). 

CORRELATION  OF  ANALYSIS  AND  EXPERIMENTS 

Ideally,  the  closed  form  solutions  for  local  delaminations  from 
matrix  cracks  should  be  verified  and  calibrated  by  comparing  them 
to  detailed  finite  element  solutions  that  sum  the  G  components  for 
the  three  fracture  modes  to  obtain  a  total  G.  Such  comparisons  have 
been  made  for  edge  delamination  [3-5].  However,  to  date  only  the 
case  of  uniform  local  delamination  from  90  degree  matrix  cracks 
has  been  verified  [8].  Part  of  the  difficulty  Is  that  the  modelling  for 


14 


local  delamination  requires  a  full  3D  analysis,  with  the  matrix  crack 
geometry  modeled  accurately.  The  solutions  for  non-90  degree 
matrix  cracks  that  have  been  attempted  to  date  [9,10],  have  assumed 
that  the  angle  ply  matrix  cracks  may  be  represented  as  straight 
cracks  through  the  laminate  thickness  that  are  normal  to  the 
interface.  However,  as  seen  in  micrographs  of  the  edges  for  the 
(O2/e2/-02)s  laminates  [1],  these  angle  ply  matrix  cracks  are  curved, 
intersecting  the  ply  interface  at  oblique  angles.  This  difference  in 
actual  versus  assumed  matrix  crack  geometry  through  the  laminate 
thickness  may  have  a  significant  effect  on  the  G  components 
calculated  from  a  finite  element  analysis  for  local  delamination. 

Although  the  fatigue  behavior  of  the  (O2/02/-02)s  laminates 
cannot  currently  be  predicted,  it  may  be  possible  to  correlate  the 
fatigue  delamination  onset  behavior  between  laminates  with 
different  values  of  0.  The  data  from  each  layup  may  be  used  to 
generate  delamination  onset  criteria  by  generating  plots  of 
maximum  cyclic  G  versus  the  number  of  cycles  to  delamination  onset 
for  the  composite  material.  A  similar  characterization  has  been  used 
previously  for  edge  delamination  data  [5,11]. 

Laminate  and  sublaminate  moduli  for  the  (O2/02/*02)s 
AS4/3501-6  graphite  epoxy  laminates  were  calculated  using  the 
properties  in  table  3  in  laminated  plate  theory.  The  moduli  of  the 
locally  delaminated  regions,  Eld.  were  calculated  using  eq(7)  to 
incorporate  the  influence  of  bending-extension  and  shear-extension 
coupling.  These  moduli  were  then  used  in  eq(2),  along  with  the 
maximum  cyclic  loads  in  the  fatigue  tests,  to  calculate  strain 
energy  release  rates  for  uniform  local  delaminations  corresponding 
to  the  measured  number  of  cycles  to  delamination  onset  for  each 
layup.  In  addition,  the  matrix  crack  length,  d,  corresponding  to  the 
number  of  cycles  for  local  delamination  onset  was  measured  from 
the  X-ray  radiographs  where  local  delamination  was  first  observed 
for  each  specimen.  Examination  of  the  partial  local  delaminations  in 
these  radiographs  indicated  that  the  geometry  shown  in  fig.  15c  was 
reasonable  for  the  layups  tested.  Therefore,  measured  matrix  crack 
lengths,  d,  were  used  in  the  expression  in  Table  4  corresponding  to 
fig.  15c  to  determine  critical  delamination  lengths  along  the  free 
edge,  a,  at  local  delamination  onset  for  each  layup.  Table  6  shows 
the  range  of  delamination  lengths  along  the  free  edge,  a,  normalized 
by  the  gage  length,  L,  measured  just  after  local  delamination  onset 
in  fatigue.  These  critical  delamination  lengths  were  used,  along  with 
the  maximum  cyclic  loads  in  the  fatigue  tests,  to  calculate  strain 
energy  release  rates  from  eq(20)  for  partial  local  delaminations 
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corresponding  for  the  measured  number  of  cycles  to  delamination 
onset  for  each  layup.  Hence,  plots  of  Gmax  versus  the  number  of 
cycles  to  delamination  onset,  N,  were  generated  for  each  layup 
assuming  either  uniform,  or  partial,  local  delamination  onset  from  a 
-0  degree  matrix  ply  crack. 

Figure  25  shows  the  Gmax  vs.  N  plot  for  both  the  uniform  and 
partial  local  delamination  models.  The  maximum  cyclic  G  levels  for 
the  uniform  and  partial  local  delamination  models  are  very  similar 
because  the  experimentally  determined  critical  delamination  lengths 
were  small  (Table  6,  fig. 20).  In  addition,  for  any  particular 
maximum  cyclic  load,  the  maximum  cyclic  G  levels  for  all  values  of 
6  also  agree  fairly  well,  indicating  that  the  assumptions  of  shear- 
extension  and  bending-extension  coupling  reflected  in  eq(7)  was 
accurate  (fig. 12).  Finally,  the  number  of  cycles  to  delamination 
onset  at  each  cyclic  G  level  appear  to  be  greater  for  the  15  and  20 
degree  cases  than  for  the  25  and  30  degree  cases.  These  differences 
in  the  number  of  cycles  to  delamination  onset  may  be  due  to  a 
difference  in  either  the  mixed-mode  ratio,  or  the  extent  of  matrix 
cracking  prior  to  delamination  onset,  for  the  individual  cases. 
Neither  of  these  effects  are  incorporated  in  the  current  analyses. 

Figure  26  shows  the  Gmax  vs.  N  plot  for  the  uniform  local 
delamination  models  with  mechanical  loading  only,  G*^,  and  with 
mechanical,  thermal,  and  hygroscopic  loading  combined,  G*^+t+h  por 
any  particular  maximum  cyclic  load,  the  maximum  cyclic  G^+t+h  is 
only  slightly  greater  than  the  the  maximum  cyclic  G^^  because  of  the 
relatively  small  contribution  of  the  residual  thermal  and  moisture 
stresses  to  G  (Table  5). 

The  Gmax  vs.  N  plots  in  figures  25  and  26  provide  some  insight 
into  the  accuracy  of  the  G  solutions  derived  in  this  paper.  However, 
as  noted  earlier,  a  detailed  analysis  of  the  G  components  and  a 
mixed-mode  fatigue  failure  criterion  for  this  material  may  be 
needed  to  predict  the  fatigue  behavior  of  these  (O2/02/-62)s 
laminates. 


CONCLUSIONS 


Constant  amplitude  tension-tension  fatigue  tests  were 
conducted  on  AS4/3501-6  graphite  epoxy  (O2/02/-92)s  laminates, 
where  0  was  15,  20,  25, and  30  degrees.  Fatigue  tests  were 
conducted  at  a  frequency  of  5  Hz  and  an  R-ratio  of  0.1.  Dye  penetrant 
enhanced  X-ray  radiography  was  used  to  document  the  onset  of 
matrix  cracking  in  the  central  -0  degree  plies,  and  the  subsequent 
onset  of  local  delaminations  in  the  0/-0  interface  at  the 
intersection  of  the  matrix  cracks  and  the  free  edge,  as  a  function  of 
the  number  of  fatigue  cycles.  Unlike  the  quasi-static  tests  [1],  the 
sequence  of  damage  was  evident  in  the  series  of  radiographs  taken 
during  the  cyclic  loading,  with  the  matrix  cracking  always  preceding 
the  onset  of  local  delamination. 

Two  strain  energy  release  rate  solutions  for  local 
delamination  from  matrix  cracks  were  derived:  one  for  a  local 
delamination  growing  from  an  angle  ply  matrix  crack  with  a  uniform 
delamination  front  across  the  laminate  width,  and  one  for  a  local 
delamination  growing  from  an  angle  ply  matrix  crack  with  a 
triangular  shaped  delamination  area  that  extended  only  partially  into 
the  laminate  width  from  the  free  edge.  Plots  of  Gmax  vs.  N  were 
generated  to  assess  the  accuracy  of  these  G  solutions.  The  influence 
of  residual  thermal  and  moisture  stresses  on  G  were  also  quantified. 
However,  a  detailed  analysis  of  the  G  components  and  a  mixed-mode 
fatigue  failure  criterion  for  this  material  may  be  needed  to  predict 
the  fatigue  behavior  of  these  laminates. 
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APPENDIX  1  -  UNIFORM  LOCAL  DELAMINATION  FROM  AN  ANGLE  PLY 

MATRIX  CRACK 

For  a  local  delamination  from  a  matrix  crack,  the  laminate 
compliance,  S,  where  S  =  1/E,  for  a  laminate  with  a  given  local 
delamination  area.  A,  in  an  interface  whose  total  area  is  A*,  may  be 
represented  as 


S-Slam  _  A 

S*-Slam  a*  (A1) 


where  S*  =  1/E*ld  .  and  E*ld  =  (tLoEioVt  as  defined  in  ref.2.  Fig.  11a 
shows  the  uniform  local  delamination  that  extends  across  the  entire 
laminate  width  and  grows  normal  to  a  90  degree  matrix  ply  crack  in 
the  lead  direction.  For  the  uniform  delamination  in  fig.  11a,  A  = 
mwa,  and  A*  =  mwL,  where  m=2  for  an  interior  matrix  crack  and  m=1 
for  a  matrix  crack  in  a  surface  ply.  Substituting  these  expressions 
into  eq(A1)  yields 


_J _ L_-)(X)a  +  -l- 

ttoELD  Elam 


(A2) 


which  is  identical  to  eq(7)  in  ref.2  for  local  delamination  from  a  90 
degree  matrix  ply  crack. 

Fig.  11b  shows  the  uniform  local  delamination  that  extends 
across  the  entire  laminate  width  and  grows  from  a  -0  degree  matrix 
ply  crack  normal  to  the  load  direction.  For  the  uniform  delamination 
in  fig.  11b,  A=  (mw/sin0)  (a  sin0)  =  mwa,  and  A*  =  mwL.  Substituting 
these  expressions  into  eq(A1)  also  yields  eq(A2). 

For  the  90  degree  matrix  crack  induced  local  delamination  in 
fig  11a,  dA  =  mwda.  Similarly,  for  the  -0  degree  matrix  crack 
induced  local  delamination  in  fig.  lib,  dA  =  (mw/sin0)  d(a  sin0)  = 
mwda.  In  ref.2,  G  is  defined  as 


2  dA 


(A3) 


Substituting  for  V  =  wLt  and  dA  =  mwda  and  differentiating  eq  (A2) 
with  respect  to  "a"  yields 
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G  =  -E?-(_J - 1 — 1 

2mw2\tLDELD  IlamElAM'  (A4) 

for  uniform  local  delaminations  from  either  a  90  degree  or  a  -6 
degree  matrix  crack. 


APPENDIX  2  -  ANALYSIS  OF  MODULI  FOR  ASYMMETRIC  SUBLAMINATES 

In  laminated  plate  theory,  the  load  and  moment  resultants  are 
related  to  the  midplane  strains  and  curvatures  as 


(A5) 


where  A,  D,  and  B  represent  the  extensional,  bending,  and  bending- 
extension  coupling  matrices,  respectively.  By  prescribing  an  applied 
mechanical  strain,  ex  =  eo,  setting  the  applied  Ny,  Nxy,  and  My  equal  to 
0,  and  setting  kx  and  Kxy  equal  to  0,  eq{A5)  becomes 


Nx 

£0 

0 

Ey 

0 

.FAB] 

Yxy 

Mx 

'  I  B  dJ 

0 

0 

Ky 

Mxy  , 

0 

(A6) 


The  three  equations  whose  left  hand  side  are  zero  may  be  written  as 

-Ai2£o  =  A22£y  +  A26Yxy  +  ^22<y 
-A  1660  =  A26ey  +  A66yxy  +  ^26^y 
-Bi2eo  =  B22ey  +  B26'yxy  +  D22Ky 


where  Ajj-Ajj.  This  set  of  three  non-homogeneous  equations  may  be 
solved  for  ey,  yxy.  and  Kxy  using  Cramer’s  rule.  Hence, 


Di  D2  Di 

ey  =  -^£o,  Y*y  =  -^eo,  Ky  =  -^Eo 


(A8) 
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where 


-A, 2 

A26 

B22 

D,= 

-Ai6 

A66 

B26 

-Bi2 

B26 

D22 

A22 

-Ai2 

B22 

D2  = 

A26 

-Ai6 

B26 

B22 

-Bi2 

t>22 

A22 

A26 

-A  12 

D3  = 

A26 

A66 

-Ai6 

B22 

B26 

-Bi2 

A22 

A26 

B22 

D  = 

A26 

A66 

B26 

B22 

B26 

D22 

These  expressions  for  ey,  Yxy.  and  xxy  may  be  substituted  into  the 
first  equation  of  eqs(A6)  to  yield 


Nx 


(A9) 


Hence,  Nx  may  be  obtained  explicitly  as  a  function  of  the  terms  in 
the  A,  B,  and  D  matrices  and  the  prescribed  axial  strain,  eq.  Then 
noting  that 


Eid  = 


2Nx 

(N-n)heo 


(A10) 


yields 


Eld  = 


-2— 

(N-n)h 


(A11) 


By  prescribing  the  loading  outlined  above,  but  setting  yxy  -  0 
instead  of  prescribing  thatNxy»  0.  i.e.,  by  prescribing  an  applied 
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mechanical  strain,  ex  »  eo,  setting  the  applied  Ny  and  My  equal  to  0, 
and  setting  yxy,  kx,  and  Kxy  equal  to  0,  eq(A5)  becomes 


Nx 

£0 

0 

Cy 

Nxy 

A  B1 

0 

Mx 

B  DJ 

0 

0 

Ky 

Mxy 

0 

(A12) 


The  two  equations  whose  left  hand  side  are  zero  are 

0  =  AiiEo  +  A22ey  +  B22Ky  (A1 3) 

0  =  BnEo  +  B22£y  +  D22Ky  (A1 4) 

where  Aij=Aji  and  Bjj=Bji.  Multiplying  eq(A13)  by  -D22,  multiplying 
eq(A14)  by  B22.  and  then  adding  the  two  equations  yields 

e  -/B12B22  -  Ai2D22tjv 

^  \A22D22  -  B22B22'  {A15) 

Furthermore,  multiplying  eq(A13)  by  *B22.  multiplying  eq{A14)  by 
A22.  and  adding  the  two  equations  yields 

^  -/A12B22  -  A22Bi2j^^ 

^  \A22D22  -  B22B22/  (A16) 


Finally,  substituting  eqs(A15)  and  {A16)  into  the  first  equation  of 
eq(A12)  yields 


Nx  = 


n” '  R  ^(a' n” ' 

VA22D22  -  D22B22'  'A22IJ22  -  D22D22' 


A22D22  •  B22B22I 


Co 


(A17) 


Hence,  Nx  may  be  obtained  explicitly  as  a  function  of  the  terms  in 
the  A,  B,  and  D  matrices  and  the  prescribed  axial  strain,  eo-  Then 
recalling  eq(A10) 


Eid 


2Nx 

(N.n)hEo 


and  substituting  the  expression  for  Nx  from  eq(A17)  yields 
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(A18) 


Eix)  = 


2 

(N-n)h 


/A12D22  -  512622)  A 
IA22D22  *  822822* 


A12B22  -  A22812 
A22D22  •  822822 


APPENDIX  3  -  MOISTURE  CONTENT  DETERMINATION 

Just  prior  to  testing  the  (02/62/-82)s  laminates,  two  control 
coupons  of  each  layup  were  weighed,  dried,  and  then  weighed  again 
to  determine  AH  for  the  laminates  in  the  ambient  environment.  These 
dried  control  specimens  were  then  weighed  periodically  to  identify 
their  moisture  content  by  weight  as  a  function  of  time. 

Fig.  27  shows  the  percentage  moisture  content  by  weight,  AH, 
as  a  function  of  the  number  of  days  from  the  time  the  specimens 
were  dried.  There  was  very  little  difference  in  the  percentage 
moisture  weight  gain  for  the  different  layups.  The  percentage 
difference  in  the  original  weight  measured  just  before  drying  and 
the  weight  measured  immediately  after  the  specimens  were  dried 
was  approximately  0.2%  for  all  the  layups.  This  value  is  shown  in 
fig. 27  as  the  original  ambient  condition.  Because  all  of  the  fatigue 
data  were  generated  within  a  period  of  30  days,  the  percentage 
moisture  weight  gain  for  all  the  laminates  that  were  tested  in 
fatigue  was  assumed  to  be  between  0.2%  and  0.3%. 
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TABLE  1 


X-ray  Schedule  for  (O2/02/-62)s  Fatigue  Tests 


Pmax.  N  (lbs) 

X-ray  taken 
every  N  cycles 

Between  Ni 
cycles 

26.690  (6000) 

1000 

1 

5000 

22,240  (5000) 

1000 

1 

5000 

17,790  (4000) 

1000 

2000 

1 

10,000 

10,000 

20,000 

15,570  (3500) 

1000 

2000 

5000 

10,000 

1 

10,000 

20,000 

50,000 

10,000 

20,000 

50.000 

110,000 

13,340  (3000) 

1  .000 

50,000 

1 

100,000 

100,000 

1,000,000 

25 


TABLE  2  -  Average  Number  of  Martrix  Cracks  at  Local 

Delamination  Onset  in  (O2/02/-62)s  AS4/3501-6 
Graphite  Epoxy  Laminates 


e,  degrees 

Pmax  .  N  (lbs) 

1  5 

20 

25 

30 

Static  Test 

1.35 

1.57 

2.67 

5.25 

17,790  (4000) 

14.2 

9.7 

11.0 

20.7 

12.0 

10.3 

13.7 

23.3 

13,340  (3000) 

9.0 

14.0 

15.7 

32.0 

TABLE  3  -  AS4/3501-6  Graphiie  Epoxy  Lamina  Properties 

Eii  =  135  GPa  (0.5  x  10^  Psi) 

E22  ®  1"*  GPa  (TB  X  10^'  Psi) 

Gi2  “5.8  GPa  (0.847  x  IO6  Psi) 

V12  =0.301 

ai  =  -0.41  X10-6/OC  (-0.23  x10-6/“F) 
a2  =  26.8  xlO-6  /“C  (14.9  x10-6  /«F) 

Pi  =  0.00  /weight  % 
p2  =  5.56  x1 0-3 /weight  % 
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TABLE  4  -  Partial  Local  Delamination  Geometries 


TABLE  5  -  Percentage  Difference  between 
and  for  AH  =  0.3% 


0,  degrees 

Pmax  .  N  (lbs) 

1  5 

20 

25 

30 

lEEEEKEIiEnHi 

5.8% 

7.5% 

7.7% 

7.0% 

15,570  (3500) 

5.0% 

6.6% 

6.8% 

6.1% 

17,790  (4000) 

4.3% 

5.8% 

6.2% 

5.8% 

22,240  (5000) 

3.7% 

4.9% 

5.1% 

4.9% 

TABLE  6  -  Range  of  Critical  Normalized  Delamination  Lengths  at  the 
Free  Edge  of  (02/92/-02)s  Graphite  Epoxy  Laminates 
Measured  After  Local  Delamination  Onset  in  Fatigue 


e.  Degrees 

a/L 

15 

.034-. 106 

20 

.036-. 090 

25 

.029-. 063 

30 

.019-. 062 

27 


Fig.  1a.  RADIOGRAPH  OF  (02/202/-202)s  LAMINATE 
AFTER  MATRIX  CRACKING  ONSET  IN  FATIGUE 
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Fig.  1b.  RADIOGRAPH  OF  (02/202/-202)s  LAMINATE 
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Matrix  crack 


S  5000  lbs 


N,  Cycles 


N,  Cycles  to  Delam 


LOCAL  DELAMINATION  MODELS 


FIG. 12  Normalized  G  for  Uniform  Local  Delamination  In 
0/  -0  Interface  of  (0/0/  -0)  Laminates 
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,  Degrees 


FIG. 13  Normalized  G  for  Edge  Delamination  in 
0/  -0  Interface  of  (0/0/  -0)  Laminates 
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,  Degrees 


Fig.  15.  PARTIAL  LOCAL  DELAMINATION  GEOMETRIES 


% 
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FIG.  16  Normalized  G  for  Delamination  in  15/-15  Interface 

of  (0/1 5/- 15)  Laminate 


Normalized  G  for  Delamination  in  20/-20  Interface 
of  (0/20/-20)  Laminate 


FIG. 18  Normalized  G  for  Delamination  in  25/-25  Interface 

of  (0/25/-25)  Laminate 


FIG. 19  Normalized  G  for  Delamination  in  30/-30  Interface 

of  (0/30/-30)  Laminate 


FIG.20  Normalized  G  for  Local  Delamination  in  0/  -0  Interface 

of  (0/0/  -0)  Laminates 
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0.02  0.04  0.06  0.08 

a/L 


FIG.21  for  local  DELAMINATION 

FROM  -1 5  DEGREE  MATRIX  CRACK 


FIG.22  for  local  DELAMINATION 

FROM  -20  DEGREE  MATRIX  CRACK 


M+T+H  pQpj  delamination 

SA  -25  DEGREE  MATRIX  CRACK 


FIG.24  for  local  DELAMINATION 

FROM  -30  DEGREE  MATRIX  CRACK 

( 0  /  3  0  /  -  3  0 )  Graphite  Epoxy 
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AS4/3501-6 


N,  Cycles  to  Delamination  Onset 
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